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Extremely high energy (~ 10 22 eV) cosmic neutrino beams initiate high energy particle cascades in 
the background of relic neutrinos from the Big Bang. We perform numerical calculations to show 
that such cascades could contribute more than 10% to the observed cosmic ray flux above 10 19 eV if 
neutrinos have ~ eV masses. The required intensity of primary neutrinos could be consistent with 
astrophysical models for their production if the maximum neutrino energy reaches to ~ 10 22 eV 
and the massive neutrino dark matter is locally clustered. Future observations of ultra high energy 
cosmic rays will lead to an indirect but practical search for neutrino dark matter. 
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It has been claimed that pure cold dark matter (CDM) 
leads to a larger baryon fraction (f It) than predicted by 
big bang nucleosynthesis (BBN) if the observed hot X- 
ray-emitting gas represents a fair sample of the universe 
0. An admixture of hot dark matter (HDM) with CDM 
shifts the estimates of the baryon fraction closer to that 
by BBN. In addition, this mixed cold + hot dark mat¬ 
ter model (CHDM) has been shown to agree well with 
the cosmic microwave background (CMB) spectrum mea¬ 
sured by COBE, and galaxy group properties such as the 
number density of clusters |Q. Neutrinos are the best 
candidate for HDM and a total neutrino mass of 5 eV, 
or m v ~ m Vr ~ 2.4 eV (fU ~ 0.2) may be a solution 
consistent with all available observations. 

If HDM consists of cosmological background neutrinos 
(CBN) with ~ eV mass and there exist cosmic neutrino 
beams reaching to ~ 10 22 eV, the interactions of these 
extremely high energy (EHE) cosmic neutrinos with the 
CBN during their propagation can become significant j|] 
due to the enhanced interaction probability at the Z bo¬ 
son resonance. The resulting neutrino cascade causes 
modifications such as a bump and a dip in the EHE neu¬ 
trino spectrum at Earth Q. The cascade contains several 
hadronic decay channels that produce mostly pions which 
in turn reproduce neutrinos through their decay |5|], but 
also 7-rays and some nucleons. Hence this mechanism 
has been proposed |||7j] as an explanation of the highest 
energy cosmic rays (HECRs) whose flux above ~ 5 x 10 19 
eV is severely attenuated by photopion production on 
the cosmic microwave background (CMB) [|| in case of 
nucleon primaries, forming the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff §, 

The evolution of cascades initiated by EHE cosmic 
neutrinos is determined by very complex chains of in¬ 
teractions: The neutrinos undergo vv reactions which 
involve the exchange of W and Z bosons and hadroniza- 


tion of their strongly interacting decay products. The 
produced photons, electrons and protons collide with the 
CMB, the infrared and optical background (IR/O) and 
the universal radio background (URB) |l(]{ll]] initiat¬ 
ing electromagnetic cascades. Electrons are also subject 
to synchrotron cooling in extragalactic magnetic fields 
(EGMF). The final particle fluxes after propagation de¬ 
pend on all these interactions and solving the relevant 
transport equations is inevitable for an accurate evalua¬ 
tion of the consequences of this scenario. In this letter 
we present numerical calculations of the “primary” EHE 
neutrinos and the “secondary” 7-rays and protons that 
may contribute a sizable fraction of the observed HECRs 
above ~ 10 19 eV, under different assumptions concern¬ 
ing neutrino mass and local density enhancement of the 
HDM. Several observable signatures to confirm or rule 
out this scenario are discussed followed by a summary. 

Cascading Calculation- Our numerical calculations 
combine simulation codes for neutrino cascades U-[l , and 
for electromagnetic cascades and nucleon propagation 
011 • Detailed accounts of these codes are provided 
in Refs. J||l2j. The following processes are included: In¬ 
elastic and elastic vv collisions involving an exchange of 
either a W or a Z boson on the CBN; the subsequent de¬ 
cays of produced 7r’s, /i’s, and r’s, hadronization of quark 
jets, all of which eventually feed into the electromagnetic, 
neutrino, and nucleon channels; 77 —> e + e _ on the CMB, 
the IR/O and the URB; inverse Compton scattering on 
the same backgrounds; triplet pair production and double 
pair production on the CMB; synchrotron cooling in the 
EGMF; the nucleon interactions on the CMB (pair pro¬ 
duction and photopion production), and neutron decay; 
redshifting and evolution of the black body temperature 
due to expansion of the universe. For the IR/O we used 
recent data & and for the (unmeasured) URB we used 
the highest prediction of Ref. [[TIJ yielding conservatively 
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low 7 -ray fluxes around 10 20 eV for which the URB is the 
most important target for pair production. We neglect 
interactions of EHE neutrinos with the CMB photons 
which are of comparable importance to those with the 
CBN only for neutrino energies above the Z resonance 



The hadronic decay of Z bosons resonantly produced 
by neutrinos of energy 

E ies = M 2 /2m„ = 4 x 10 21 (^) 1 eV, (1) 

with the CBN is the most important neutrino process for 
production of y-rays and nucleons whose spectra are de¬ 
termined by the hadron fragmentation function. At the 
energy range around the Z pole, this has been measured 
accurately by the LEP at CERN. We implemented em¬ 
pirical functions using the MLLA approximation @ in 
our code, which have been fitted by measurements of the 
inclusive production rates of tt 31 and pp with the OPAL 
detector 0 - This constitutes the major revision of the 
original codes in Refs. ||||. 

The dominant contribution to secondary particle fluxes 
from resonant Z production can be estimated analyti¬ 
cally, for example for the produced y-ray spectrum, 
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where Tz — 0.03 Mz is the decay width of the Z 
boson whose mass is Mz , drih/dx is the 7 r° frag¬ 
mentation spectrum, and A z — 38 Gpc is the mean 
free path of neutrinos at E Tes given by Eq. (©■ 
The number of photons above 3 x 10 19 / (m v / eV) eV 
is then 0.012 E les (dN u /dE u )(E res ) Mpc -1 , compared to 
0.03 from our full numerical calculation which includes 
contributions from all channels and uses the more accu¬ 
rate MLLA formula for the fragmentation function. 

The particle fluxes - For a general discussion we con¬ 
sider a homogeneous distribution of sources radiating 
EHE neutrinos with a constant differential spectrum 
oc E~ q and a luminosity per comoving volume that scales 
as (1 + z) m between z = z m - ln and z = z max (dj, with 
m characterizing source evolution. Because of the small 
neutrino absorption probability, the results are essen¬ 
tially independent of z m i n < 1. We assume a flat uni¬ 
verse with a Hubble constant of Hq = 65 km sec -1 
Mpc -1 which is consistent with the CHDM picture of 
the universe [jj]j2j]. We use q = 1, a typical spectral index 
expected for neutrinos produced from photopion produc¬ 
tion by accelerated protons (ll] | 2 (]]. Results are, however, 
rather insensitive to q for q < 2. The ratio of emitted 
and v e fluxes is assumed to be ~ 1 . 86 , as expected from 
charged pion decay. Furthermore, for the moment we as¬ 
sume that the source luminosity in y-rays and nucleons 
is negligible compared to the neutrino luminosity. 

HDM is usually expected to cluster locally and 
the Fermi distribution with a velocity dispersion v 


yields an overdensity f v < v 3 m 3 /(2Tt) 3 l 2 jn v ~ 

330 (t>/500kmsec -1 ) 3 (m v /eV) 3 over the uniform CBN 
[fH| . If clustering occurs on a scale /„ smaller than the 
typical attenuation length Z att ~ few Mpc of nucleons and 
y-rays around 10 2 ° eV, the ratio of their fluxes produced 
on that scale to the ones produced on the uniform back¬ 
ground is ~ fJv/htf Therefore, while clustering in the 
galactic halo or in a nearby galaxy cluster is unlikely to 
contribute to the HECR flux ^2| , neutrinos clustering in 
the local Supercluster may have /„ ~ 100, l„ ~ few Mpc. 
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FIG. 1. Energy spectra of nucleons, y-rays and neutrinos 
for the scenario described in the text. Error bars are the 
combined data from the Fly’s Eye |2(| and the AGASA plj ] 
experiments above 10 19 eV. Also shown are piecewise power 
law fits to the observed charged CR flux below 10 19 eV, the 
EGRET measurement of the diffuse y-ray flux between 30 
MeV and 100 GeV, and experimental neutrino flux limits from 
Frejus and Fly’s Eye [^, as well as projected neutrino 
sensitivities of the future Pierre Auger and NASA’s OWL 
[p9| projects. 


In Fig. [j] we show the calculated spectra for the fol¬ 
lowing typical case: m = 3, z m i n = 0, z max = 3, an 
EGMF < 10 -12 Gauss, m Ve = m VlM = m Vr = 1 eV, 
f v — 300, and l v = 5 Mpc. It can be seen clearly 
that the predicted fluxes are consistent with the mea¬ 
surement of the diffuse y-ray flux by EGRET [23] and 
with upper limits on neutrino fluxes by Frejus 24] and 
Fly’s Eye ]25|]. Typically, the energy content in the pro¬ 
duced low energy cascade y-rays is a few percent of the 
neutrino energy which agrees with a rough analytical es¬ 
timate giving ~ 10/(H o \z)Tz/Mz- By scaling the cos- 
mologically produced low energy y-ray flux in Fig. [j] 
with la.tt/(fvlv), the EGRET constraint on diffuse y-rays 
requires /„ > 20 (G/5Mpc) -1 . 

The EHE part of the secondary y-rays and pro¬ 
tons possibly constitute a hard component of the ob¬ 
served HECRs without a GZK cutoff. The energy 
content in this “visible” HECR component is about 
(T z /M z )(fvlv/\z)E/ es (dNv/dEv)(E les ), again consis¬ 
tent with the fluxes shown in Fig. |l]. The collisions of the 
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EHE cosmic neutrinos with the HDM can be responsible 
for ~ 10% of the observed cosmic rays above 10 19 eV with 
dominant contributions of 7 -rays above the GZK cutoff. 
The fluxes deviate at most by 50% for m„ e <C . 



FIG. 2. Same as in Fig. [I] but for the case of an overdensity 
of 20 over 5 Mpc, showing fluxes above 10 18 eV. Line key as 
in Fig. |l] with the “visible” sum of 7 -ray and nucleon fluxes 
shown as thin solid line in addition, as well as upper limits in 
the bins where no HECR were seen below the highest energy 
event. This case sets an upper bound for the intensity of 
primary EHE neutrino beams allowed by the EGRET diffuse 
7 -ray limit, assuming only secondaries of neutrino interactions 
contribute to the EGRET flux. 

Fig. H shows the high energy part of the resultant 
spectra above 10 18 eV as in Fig. [lL but for the case of 
the lower local enhancement of the neutrino dark matter, 
/„ = 20 over a scale l v = 5 Mpc, the lowest possible f v 
allowed by the EGRET bound. As compared to the case 
for stronger clustering shown in Fig. [l], the required EHE 
neutrino intensity is 10 times larger. 
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FIG. 3. Same as in Fig. ^ but for the optimistic case of 
an overdensity of 10 3 over 5 Mpc that would be required if 
L 7 ~ 13Li,/3. 


In general, models based on photopion production pre¬ 
dict an integrated photon source luminosity L 1 that 
is comparable to the total neutrino luminosity, L 7 ~ 
|li],|3(|. In this case, the EGRET constraint 
translates into the more stringent requirement f v > 
10 3 (Z„/5Mpc) _1 , as can be seen by applying the above 
mentioned scaling to the integrated neutrino luminosity 
from Fig. ^j. Fig. || shows the fluxes for this optimistic 
case of strong clustering. This bound on f v can be re¬ 
laxed if most of L-y does not appear at EGRET energies, 
but is dominantly released in the TeV range. This could 
be a detectable signature from individual point sources 
& in addition to the secondary 7 -rays from neutrino in¬ 
teractions appearing at EGRET energies. Furthermore, 
the scenario discussed here requires sources that are opti¬ 
cally thick for accelerated protons with respect to photo¬ 
pion production because otherwise the observable proton 
flux below the GZK cutoff would be comparable to the 
neutrino flux | 2 (]]. 

Discussion The EHE neutrino scenario we explored 
here is quite solid in terms of the particle physics because 
the interactions with the cosmological backgrounds oc¬ 
cur in the well measured LEP energy range. No physics 
beyond the Standard Model is involved except neutrino 
mass. The major uncertainty arises in the question 
whether any astronomical objects are capable of produc¬ 
ing neutrinos with energies of ~ E les . In the conven¬ 
tional models invoking the decay of photoproduced pi- 
ons, primary protons must be accelerated to ~ 20E res ~ 
10 23 /(m„/eV) eV in order to generate neutrinos of en¬ 
ergy Fires- Furthermore, the sources would need a dense 
photon target to supply high neutrino luminosity and to 
absorb protons and 7 -rays. Thus, a new model for the 
neutrino beam sources may be necessary |3^]. 

Interestingly, the energy generation rate of ~ E les 
neutrinos for the scenario shown in Fig. |lj. 1.8 x 10 45 
erg Mpc ~ 3 yr _1 divided by the the rate of cosmological 
Gamma Ray Bursts (GRBs), 3 x 10 -8 Mpc -3 yr _1 [ fl3| , 
yields ~ 6 x 10 52 erg, and is comparable to the observed 
energy release including afterglow from a typical GRB in 
the BATSE range [|dj. 

The EHE neutrino scenario has several advantages to 
explain the HECR observation. The observed relatively 
hard spectrum without GZK cutoff 271 is reasonably re¬ 
constructed in our model which is determined mainly by 
the well-measured hadron fragmentation function at the 
Z pole and the energy loss process in the cosmological 
backgrounds, regardless of the nature of the EHE neu¬ 
trino sources. The highest energy events above the GZK 
cutoff can originate from very distant powerful objects 
because neutrinos propagate without significant energy 
loss. For example, the AGN 3C147 at a redshift of 0.545 
is a candidate for the Fly’s Eye 3x 10 2 ° eV event |^5|. For 
the same reason, it is natural that we found no nearby 
powerful astronomical objects in directions of the pos¬ 
sible event clusters observed by AGASA 1361. Because 
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the EHE neutrino beams can be responsible for a sizable 
fraction of cosmic rays above 10 19 eV, this scenario can 
explain the observational fact that the intensity of the 
events observed above the GZK cutoff is consistent with 
the extrapolation of the flux from lower energies. 

Among the observable signatures of the neutrino sce¬ 
nario are the primary EHE neutrinos whose flux should 
be detectable, as projected sensitivities of future exper¬ 
iments such as the Pierre Auger ^ 8 | and NASA’s OWL 
p9f projects suggest. Correlation of the arrival direction 
of the EHE neutrino and the secondary HECR showers 
may also be observable. As opposed to conventional mod¬ 
els with nucleon primaries, our model predicts that some 
of the observed HECRs should originate in sources at 
cosmological distances. At energies beyond the GZK cut¬ 
off, the correlation of arrival directions of HECR showers 
with sources at distances t should be easy to de¬ 
tect since the background from a conventional nucleon 
component should be suppressed due to the GZK effect, 
whereas the component proposed here continues as a rel¬ 
atively flat spectrum. Finally, this scenario predicts a 
7 -ray domination above 10 2 ° eV, and next generation 
experiments should settle the question whether observed 
HECR are consistent with 7 -ray primaries [ fpif . 

In summary, we have seen that collisions of EHE cos¬ 
mic neutrino beams with ~ eV mass neutrino dark mat¬ 
ter would explain the observed HECR energy spectrum, 
regardless of the nature of the neutrino sources if the 
maximum neutrino energy reaches to the Z boson pole 
region and the dark matter is clustered on the Super¬ 
cluster scale by amounts consistent with expectations. 
Although EHE neutrino sources require very high effi¬ 
ciency of converting the energy to neutrino flux which 
may require a new production mechanism of neutrinos, 
the necessary neutrino intensity can be consistent with 
observed diffuse 7 -ray fluxes and the CRB energy release 
rate. The EHE neutrino scenario is a way of producing a 
relatively flat component of nucleons and 7 -rays that pro¬ 
vides a significant fraction of the HECR flux above 10 19 
eV, dominating above the GZK cutoff, without invok¬ 
ing physics beyond the Standard Model except neutrino 
mass. Future observations of HECRs lead to indirect 
search for signatures of neutrino dark matter. 
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